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ABSTRACT 

Large  quantities  of  explosives,  frequently  exceeding  1 00,000  lb,  in  different  types  of 
munitions  may  be  stored  in  a  single  magazine.  The  storage  of  munitions  has  always 
presented  safety  problems,  and  various  regulations  have  been  developed  over  the  years  to 
ensure  safe  storage  practices.  In  an  effort  to  increase  magazine  limits,  it  has  been  proposed 
that  the  maximum  credible  event  in  an  accident  scenario  may  be  significantly  reduced  if  the 
munition  store  is  divided  into  two  or  more  stacks  of  ammunition  separated  by  barriers 
designed  to  prevent  propagation  of  an  explosion  from  one  stack  to  another.  A  combined 
analytical  and  experimental  study  was  proposed  to  assess  this  hazard  and  to  determine 
whether  such  barriers  can  be  designed.  The  simulations  were  used  to  determine  the  velocity 
of  a  sand  barrier  on  impact  with  an  acceptor  ammunition  stack  in  order  to  design  meaningful 
experiments.  The  donor  ammunition  stack  was  simulated  by  a  volume  of  bare  explosive.  The 
initial  position  and  size  of  the  donor  charge  were  varied  in  different  computations.  The 
thickness  of  the  barrier  was  also  varied.  Velocity,  pressure,  and  impulse  histories  were 
monitored  at  several  stations  in  and  near  the  barrier.  Computational  results  show  that  the 
kinetic  energy  imparted  to  a  barrier  decreases  with  its  thickness,  indicating  that  thin  fast- 
moving  barriers  have  a  potential  to  do  greater  damage  to  ammunition  in  an  acceptor  stack 
than  thick  siow-moving  barriers.  Thus,  the  barriers  must  be  designed  thick  enough  to  prevent 
fragment  penetration. 


183 


Report  Documentation  Page 


Form  Approved 
0MB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  0MB  control  number. 


1.  REPORT  DATE 

AUG  1992 


2.  REPORT  TYPE 


3.  DATES  COVERED 

00-00-1992  to  00-00-1992 


4.  TITLE  AND  SUBTITLE 

Two-Dimensional  Simulations  of  Sand  Barrier  Motion  Induced  by  The 
Explosion  of  an  Ammunition  Stack  Inside 

6.  AUTHOR(S) 


7.  PEREORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  Ballistic  Research  Laboratory ,LABCOM, Aberdeen  Proving 
Ground, MD, 21005-5066 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 


11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

See  also  ADA260984,  Volume  I.  Minutes  of  the  Twenty-Fifth  Explosives  Safety  Seminar  Held  in  Anaheim, 
CA  on  18-20  August  1992. 


14.  ABSTRACT 

see  report 


15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 


a.  REPORT 

unclassified 


b.  ABSTRACT 

unclassified 


c.  THIS  PAGE 

unclassified 


17.  LIMITATION  OF 

18.  NUMBER 

ABSTRACT 

OE  PAGES 

Same  as 

30 

Report  (SAR) 

19a.  NAME  OE 
RESPONSIBLE  PERSON 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


ACKNOWLEDGMENTS 


This  work  was  supported  by  tiie  Project  Manager  for  Ammunition  Logistics,  Picatinny 
Arsenal,  NJ.  Mr.  O.  Lyman  of  the  U.S.  Army  Ballistic  Research  Laboratory,  who  was 
responsibie  for  the  experimental  study,  was  instrumental  in  providing  guidance  on  the 
magazine  design  and  munitions  stack  arrangements  which  were  used  in  this  study. 


184 


1.  INTRODUCTION 


Large  quantities  of  explosives,  frequently  exceeding  1 00,000  lb,  in  different  types  of 
munitions  may  be  stored  in  a  single  magazine.  The  storage  of  ammunition  has  always 
presented  safety  problems,  and  various  regulations  have  been  developed  over  the  years  to 
ensure  safe  storage  practices  (DOD  6055.9  STD  1984;  AR  385-64  1987).  The  situation  is 
further  complicated  in  areas  where  increased  population  density  and  reassessment  of  hazards 
have  reduced  the  applicable  explosive  limits  for  existing  magazines  (Lyman  1983).  To 
reduce  the  total  land  requirements  for  magazine  facilities,  it  is  desirable  to  limit  the  size  of  the 
maximum  credible  event  to  about  50,000  lb  (further  reductions  don’t  reduce  the  safe  inhabited 
building  distance). 

It  has  been  suggested  that  the  maximum  credible  event  in  an  accident  scenario  may  be 
significantly  reduced  if  the  munitions  store  is  divided  into  two  or  more  stacks  of  ammunition 
separated  by  barriers  designed  to  prevent  propagation  of  an  explosion  from  one  stack  to 
another  caused  by  fragments  from  the  donor  stack.  For  typical  stacks  of  ammunition,  this 
may  require  a  wall  of  sandbags  about  1  m  thick.  When  using  such  a  barrier  to  intercept 
fragments,  the  problem  of  the  effect  of  its  impact  on  the  acceptor  stack  arises.  While  the 
fragment  hazard  may  be  eliminated,  the  moving  bam'er  may  cause  sufficient  damage  to  the 
munitions  in  the  acceptor  stack  to  initiate  the  explosive  they  contain,  thus  defeating  its 
intended  purpose. 

A  combined  analytical  and  experimental  investigation  was  proposed  to  assess  this  hazard 
and  to  determine  whether  such  barriers  can  successfully  be  designed.  The  analytical  study, 
which  has  been  completed,  included  three  numerical  simulations  using  the  HULL  code.  The 
first  simulation  was  used  to  determine  the  velocity  of  the  barrier  on  impact  with  the  acceptor 
stack  in  the  magazine  environment  in  order  to  design  meaningful  experiments.  In  addition, 
several  barrier  design  issues  were  addressed  analytically.  The  second  simulation  was  used  to 
determine  the  experimental  configuration  required  to  produce  the  desired  barrier  velocity.  The 
final  simulation  was  used  to  assess  the  loading  on  typical  acceptor  ammunition  that  might  be 
antic  ipated,  in  order  to  determine  instrumentation  requirements. 
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2.  DESCRIPTION  OF  THE  MAGAZINE  AND  AMMUNITION  STACKS 


Munitions  must  be  stored  such  that  they  can  be  easily  accessed,  inspected,  and 
inventoried.  Dissimilar  items  are  stored  together,  and  each  type  must  be  available  for  removal 
without  disturbing  other  stored  munitions.  Much  of  the  munitions  stockpile  is  stored  in 
standard  magazines.  Elevation  and  plan  views  of  such  a  structure  including  a  proposed  stack 
arrangement  are  shown  in  Figure  1 . 

The  magazine  is  constructed  with  reinforced  concrete  walls,  roof,  and  floor  and  is  covered 
with  earth.  The  interior  of  the  magazine  Is  semicyiindrical,  24.38  m  long  and  8.00  m  in 
diameter.  A  steel  door,  1 .22  m  (4  ft)  wide  by  2.44  m  (8  ft)  high,  is  located  in  the  middle  of  one 
of  the  end  walls.  The  door  is  used  for  access  in  order  to  store  or  remove  ammunition  and 
also  serves  as  a  vent  in  case  of  an  explosion  inside  tiie  magazine. 

Using  palletized  M107  155-mm  projectiles,  it  is  possible  to  store  506  pallets  in  a  stack 
3.35  m  deep,  reaching  to  the  roof  of  the  magazine.  The  high  explosive  weight  per  pallet  is 
approximately  55.9  kg  (123  lb),  which  gives  28,290  kg  (62,238  lb)  per  stack.  Up  to  four  such 
stacks  could  be  stored  in  a  magazine,  while  still  leaving  room  for  barriers  up  to  1  m  thick  and 
1  m  of  access  space  on  each  side  of  each  stack. 

3.  DESCRIPTION  OF  THE  PROPOSED  EXPERIMENT 

It  Is  clearly  impractical  to  conduct  full-scale  tests  of  candidate  magazine  configurations  in 
order  to  find  an  arrangement  which  produces  the  desired  results.  As  an  alternative,  tests  in 
which  a  barrier  segment  is  launched  at  a  representative  velocity  toward  an  acceptor  stack 
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Figure  1 . 


End  Elevation  and  Plan  Views  of  a  Typical  Magazine  Showing  a  Proposed  Stacking 
Arrangement. 
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segment  have  been  proposed,  in  these  tests,  the  barrier  thickness  remains  fuii-scaie,  but  its 
laterai  extent  is  iimited. 

The  test  setup  is  iiiustrated  schematicaiiy  in  Figure  2.  in  order  to  account  for  the  iaterai 
confinement  provided  by  the  fuii  magazine  environment,  the  tests  are  conducted  in  a  trench 
2  m  deep,  3  m  wide,  and  about  10  m  iong  with  a  2-m-thick  earth  overburden.  The 
configuration  is  symmetricai,  aiiowing  two  barriers  to  be  launched  toward  two  acceptor  stacks 
by  a  centrai  ammonium  nitrate  -  fuei  oii  (ANFO)  charge.  The  acceptor  stacks  inciude  both  iive 
and  inert  155-mm  ammunition. 

4.  BRiEF  DESCRiPTiON  OF  THE  HULL  CODE 

We  used  the  HULL  code  to  conduct  the  an^yticai  study.  The  HULL  system  consists  of 
programs  for  generating  and  soiving  two-  and  tiiree-dimensional  dynamic  continuum 
mechanics  probiems  in  Euierian  and/or  Lagrangian  frameworks  pius  many  peripherai 
programs  required  to  provide  graphicai  renditions  of  the  resuits.  Numerous  constitutive 
modeis  are  inciuded,  and  materiai  parameters  are  provided  in  the  HULL  iibrary.  Detonation  is 
treated  by  the  programmed  bum  method.  We  used  the  two-dimensionai  Euierian 
finite-difference  capabiiity  (without  heat  conduction  or  viscosity)  for  our  simuiation  and  the 
graphics  routines  to  obtain  pressure  and  density  contour  piots  as  weii  as  history  piots  (Huii 
User’s  Ciass). 

5.  DESCRiPTiON  OF  THE  MAGAZiNE  SiMULATiON 

in  order  to  keep  the  size  of  the  simuiation  manageabie,  we  used  a  two-dimensionai 
axisymmetric  representation  of  the  magazine  as  iiiustrated  in  Figure  3,  The  simuiated 
magazine  is  thus  a  compiete  cyiindrical  shell  structure  having  twice  the  interior  volume  of  the 
actual  semicylindrical  magazine.  Because  of  the  absence  of  any  angular  motion  component, 
a  plane  of  symmetry  exists  at  tee  magazine  floor  (or  any  similar  plane),  effectively  providing  a 
perfectly  rigid  surface. 
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Figure  2.  Side  Elevation  and  Plan  Views  of  the  Proposed  Experimental  Arrangement. 
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ALL  MEASUREMENTS:  tn 


SAND 


A  concrete  shell  and  end  walls  enclose  the  air-filled  interior  of  the  magazine,  and  a  sand 
layer  simulates  the  earth  overburden.  Air  also  fills  the  region  surrounding  the  magazine.  No 
attempt  was  made  to  represent  the  doors. 

The  donor  ammunition  stack  is  simulated  by  a  cylinder  of  bare  TNT  representing  (in  the 
baseline  case)  56.5  Mg  (about  124,000  lb)  or  twice  the  actual  mass  (corresponding  to  the 
doubling  of  the  interior  volume).  The  initial  position  of  the  donor  charge  and  its  mass  was 
varied  in  different  computations.  In  each  case,  it  was  initiated  on  the  axis  at  the  center  of  the 
charge. 

A  sand  barrier  is  placed  near  the  center  of  the  magazine.  Sand  wall  thicknesses  varying 
from  1  m  to  3  m  were  modeled  in  different  computations.  In  a  few  computations,  vents  (or 
open  spaces)  were  left  in  the  barriers  to  facilitate  pressure  equilibration  across  them.  The 
vents  were  placed  1 .0,  2.2,  and  3.4  m  above  the  floor.  The  size  of  the  vents  ranged  from 
0.20  to  0.60  m.  Several  stations  in  and  near  the  sand  barrier  were  chosen  for  monitoring 
velocity,  pressure,  and  impulse  histories.  These  are  shown  in  Figure  4. 

The  computational  region  was  discretized  into  350  zones  covering  28  m  in  the  axial 
direction  and  75  zones  covering  6  m  in  the  radial  direction. 

Equations  of  state  for  air,  concrete,  sand,  and  TNT  products  (Jones-Wilkins-Lee  [JWL]) 
were  taken  from  the  materials  library  of  the  HULL  code. 

6.  RESULTS  OF  THE  MAGAZINE  SIMULATION 

Figures  5a  through  5e  shows  a  series  of  pressure  and  density  contour  plots  for  a  typical 
computation.  Pressure  contours  are  shown  on  the  left  side  of  the  axis  and  density  on  the 
right.  After  the  TNT  detonated,  it  took  about  3  to  4  ms  for  the  blast  wave  to  arrive  at  the  sand 
barrier.  Acceleration  of  the  wall  took  place  over  a  considerably  longer  period  of  time.  Each 
computation  was  allowed  to  run  for  30  to  40  ms.  This  was  more  than  enough  time  to  produce 
a  steady  wall  velocity.  The  contour  plots  show  increasing  diffusion  of  the  material  boundaries 
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Figure  5a.  Pressure'  and  Density  Contour  Plots  for  a  Representative 
Magazine  Simulation  Computation.  0  ms. 
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Figure  5c.  Pressure  and  Density  Contour  Plots  for  a  Representative 
Magazine  Simulation  Computation.  4  ms. 
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Figure  5e.  Pressure  and  Density  Contour  Plots  for  a  Representative 
Magazine  Simulation  Computation.  20.0  ms. 
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at  late  time.  Corresponding  plots  of  velocity  history  at  several  of  the  stations  are  shown  in 
Figure  6.  Somewhat  different  steady-state  values  are  produced  at  the  different  stations.  A 
total  of  18  computations  were  completed.  The  steady-state  velocities  computed  at  stations  7, 
8,  and  9  (near  the  leading  edge  of  the  barrier)  are  summarized  in  Table  1  along  with  the 
average  velocity  at  all  the  stations  in  the  barrier.  Some  variation  in  the  velocity  from  station  to 
station  is  evident,  but  the  general  trends  remain  dear. 

As  the  barrier  thickness  is  increased  while  retaining  a  separation  of  5.0  m  between  the 
center  of  the  charge  and  the  center  of  the  barrier  with  a  28.2-Mg  charge,  the  resulting  terminal 
velocity  decreases.  Average  terminal  velocity  is  plotted  as  a  function  of  barrier  thickness  in 
Figure  7a.  As  the  separation  distance  is  increased,  the  terminal  velocity  decreases  as  shown 
in  Figure  7b.  Approximate  representations  of  the  bam'er  momentum  and  kinetic  energy  per 
unit  area  may  be  obtained  by  multiplying  the  initial  barrier  mass  per  unit  area  by  the  velocity 
and  half  the  square  of  the  velocity,  respectively.  These  are  plotted  as  functions  of  barrier 
thickness  in  Figure  8.  The  results  are  most  consistent  with  constant  momentum  and  show 
decreasing  kinetic  energy  with  increasing  thickness. 

Venting  of  the  barriers  left  their  terminal  velocities  virtually  unchanged  as  is  evident  in 
Table  1. 

The  distances  between  the  barrier  and  the  donor  diarge  and  the  mass  of  the  donor 
charge  were  also  varied  in  different  computations.  Blast  scaling  laws  suggest  that  the  impulse 
delivered,  and,  hence,  the  terminal  momentum  of  the  barrier,  is  a  function  of  the  scale 
distance  defined  by  dividing  the  physical  distance  by  the  cube  root  of  ttie  charge  mass.  The 
variations  in  approximate  terminal  momentum  as  a  function  of  scale  distance  are  shown  in 
Figure  9.  The  results  indicate  that  this  relationship  is  reasonably  well  followed  especially  at 
larger  scale  distances. 
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Table  1.  Summary  of  Magazine  Computations 


Charge 

Mass 

(Mg) 

Barrier 

Thickness 

(m) 

Center-to-Center 

Separation 

(m) 

Terminal  Velocity 
at  Station 

7  8  9 

(m/s)  (m/s)  (m/s) 

Average 

Terminal 

Velocity 

(m/s) 

19.5 

2.0 

5.0 

■■■■ 

200 

170 

208 

28.2 

0.5 

5.0 

800 

867 

28.2 

1.0 

5.0 

500 

460 

514 

28.2 

1.0 

10.0 

440 

385 

510 

437 

28.2 

1.0“ 

10.0 

420 

400 
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443 

28.2 

1.0 

11.0 

450 

400 

530 

415 

28.2 

2.0 

5.0 

320 

290 

255 

290 

28.2 

2.0 
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300 

280 

235 

282 
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2.0 

7.5 

225 

230 

190 

224 
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10.5 

240 

240 

230 

210 
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2.0“ 

10.5 

220 

240 
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150 

150 
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140 
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142 
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Figure  7a.  Average  Terminal  Velocity  as  a  Function  of  Barrier  Thickness  for  a 
Charge  Center  to  Barrier  Center  Separation  Distance  of  5.0  m. 
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Figure  7b.  Average  Terminal  Velocity  as  a  Function  of  Separation  Distance. 
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of  Barrier  Thickness  for  a  Charge  Center  to  Barrier  Center  Separation 
Distance  of  5.0  m. 
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Figure  9.  Approximate  Terminal  Momentum  as  a  Function  of  Scale  for  All  Computation. 
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7.  DESCRIPTION  AND  RESULTS  OF  THE  TEST  SIMULATION 


A  simulation  of  the  experimental  configuration,  as  illustrated  in  Figure  1 0,  was  also  made. 
Axisymmetry  was  deemed  the  best  way  to  represent  the  lateral  confinement  in  a 
two-dimensional  simulation.  Additional  symmetry  is  afforded  by  the  experimental  setup  in 
which  identical  sand-filled  barriers  are  arranged  on  either  side  of  a  TNT  or  ANFO  charge 
which  fills  the  space  between  them.  The  computational  region  thus  extends  from  the  center  of 
the  charge  (where  a  reflective  plane  is  placed)  to  a  point  2  m  beyond  the  end  of  the  trench. 
The  thickness  of  the  barrier  was  maintained  at  1  m  in  all  the  computations,  but  the  thickness 
of  the  charge  was  varied  from  0.30  to  0.40  m  in  different  computations.  The  velocity  of  the 
wall  was  monitored  at  three  locations  (r=0.50,  1 .00,  and  1 .50  m)  along  the  radial  axis,  as 
shown  in  Figure  10. 

Representative  velocity  histories  are  shown  in  Figure  1 1 .  Because  the  driving  explosive  is 
initially  in  contact  with  the  barrier,  the  barrier  does  not  quite  reach  terminal  velocity  before 
striking  the  end  of  the  trench.  The  maximum  velocities  determined  in  four  computations  are 
summarized  in  Table  2. 
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ALL  MEASUREMENTS:  m 


Figure  10.  Experiment  Simulation  Configuration  Showing  the  Locations  of  the 
Laaranoian  Stations. 
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Figure  11.  Velocity  Histories  at  Several  Laaranqian  Stations  for  a  Representative  Test 
Simulation  Computation. 
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Table  2.  Summary  of  Test  Simulation  Computations 


Charge 

Explosive 

Thickness 

(m) 

Maximum 
Velocity  at  Station 

1  2  3 

(m/s)  (m/s)  (m/s) 

Average 

Velocity 

(m/s) 

TNT 

0.20 

270 

280 

320 

290 

TNT 

0.30 

450 

450 

480 

460 

ANFO 

0.30 

320 

320 

350 

330 

ANFO 

0.40 

400 

400 

450 

417 

8.  DESCRIPTiON  AND  RESULTS  OF  THE  ACCEPTOR  RESPONSE  SiMULATION 

A  final  calculation  was  made  to  determine  the  pressure  that  might  be  expected  inside  a 
155-mm  M107  artiliery  round  subjected  to  barrier  impact  at  360  m/s.  The  problem 
configuration  inciuding  the  iocations  of  the  pressure  monitoring  stations  is  shown  in  Figure  12. 
The  computed  pressure  histories  are  shown  in  Figure  13.  The  pressures  observed  were 
generally  lower  than  0.35  GPa. 

9.  SUMMARY  AND  CONCLUSIONS 

The  fact  that  the  kinetic  energy  imparted  to  a  barrier  decreases  with  its  thickness  indicates 
that  thin  fast-moving  barriers  have  the  potential  to  do  greater  damage  to  ammunition  in  an 
acceptor  stack  than  thick  slow-moving  barriers.  Barriers  must  be  designed  thick  enough  to 
prevent  fragment  penetration.  Additionai  thickness,  where  possible,  is  also  desirable. 

The  terminal  momentum  of  the  barriers  was  found  to  be  roughly  constant  for  a  given 
charge  mass  and  iocation.  The  application  of  biast  scaling  laws  showed  that  these  might  be 
useful  as  predictive  tools  when  ttie  distance  between  stack  and  bamer  is  not  too  small. 
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Figure  12.  Munition  Response  Simulation  Configuration  Showing  the  Locations  of  the  Laqranqian  Stations. 


Figure  13.  Pressure  Histories  at  the  Laqranqian  Stations  for  the  Munition  Response 
Simulation. 
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Terminal  velocities  were  not  significantly  reduced  when  vents  were  placed  in  the  barriers. 
Thus,  these  computations  provide  little  support  for  this  technique  as  a  method  for  mitigating 
the  hazard  posed  by  moving  barriers. 
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